Spin-Orbit Effects in Diluted Magnetic Semiconductors 



Byounghak Lee 1 , Xavier Cartoixa 1 , Nandini Trivedi 2 , and Richard M. Martin 1 
1 Loomis Laboratory of Physics, University of Illinois at Urbana- Champaign, 
1110 W. Green Street, Urbana, Illinois 61801 and 
2 Department of Theoretical Physics, Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400005 INDIA 

(Dated: February 2, 2008) 

We present a theoretical study of diluted magnetic semiconductors that treats the local sp-d 
exchange interaction J between the itinerant carriers and the Mn d electrons within a realistic band 
structure and goes beyond previous mean-field approaches. In case of Gai-^Mn^ As, we find that 
strong exchange potentials tightly bind carriers near impurity sites. Spin-orbit coupling is found 
to have a more pronounced effect on spin polarization at weak coupling and this feature can be 
exploited for magnetotransport. For low doping regime, we predict that spin-orbit coupling splits 
impurity bands and can induce a novel insulating phase. 
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It has already been more than half a decade since 
ferromagnetism with Curie temperatures T c ~ 100 K 
was discovered in III-V diluted magnetic semiconduc- 
tors (DMS) 0, and following that a large amount 
of experimental knowledge has been accumulated [2|. 
There also has been a g reat deal of theoretical activ- 
ity H II H H II H H Efl EH to understand the mecha- 
nism of ferromagnetism, but consensus has not yet been 
reached. The generally accepted theoretical picture is 
that the long range ferromagnetic order between the mag- 
netic impurities with large magnetic moment is mediated 
via the interaction with the itinerant carriers |l2i |. The 
comprehensive theoretical understanding, however, has 
been hampered by the complexity of the system arising 
from the interplay of equally important but distinct ef- 
fects, namely the electron-local moment interaction, dis- 
order effects, and spin-orbit (SO) effects. 

Transport properties indicate that DMS do not 
belong to either of the well established classes of mag- 
netism, i.e., itinerant ferromagnetism and local moment 
ferromagnetism, but rather are intermediate in charac- 
ter. Disorder plays the major role in determining the 
nature of the metal- insulator transition 0, In a d- 
dition, in typical DMS, the orbital degrees of freedom 
couple to the spin degrees of freedom via relativistic ef- 
fects. In III-V semiconductors, the SO coupling ranges 
from a few tens to several hundreds of meV, in particular 
the splitting of the bands at k=0 is 0.34 cV for GaAs 
This energy is as big as the Fermi energy in the relevant 
impurity doping concentrations. The approaches based 
on itinerant ferromagnetism often assume a continuum 
distribution of magnetic impurities interacting with host 
valence band electrons 0, 13, 13 • While this gives good 
descriptions in the metallic regime, it is not applicable 
when the local character of the magnetic interaction is 
important, not to mention the disorder effects that are 
ignored completely. More realistic tight-binding models 
based on atomic orbitals have been limited to the dilute 
impurity regime |ll| . The alternative approach of lo- 
cal moment magnetism has neglected spin-orbit coupling 
and used only a single band for each spin species @,H,Ui 



largely due to the computational complexity. This, how- 
ever, cannot be justified without an obvious symmetry 
breaking which completely lifts the orbital degeneracy. 
Density functional theory (DFT) calculations lUISI suf- 
fer high computational cost in dealing with systems with 
dilute impurities, and more complicated issues such as 
disorder effects are hard to study using small systems. 
The purpose of this paper is to introduce a robust model 
which can be applied to wide ranges of material param- 
eters without ignoring any of the important effects men- 
tioned above, and, as an immediate example, to explore 
the SO effects on the ground state electronic structures. 

Our main results are as follows: (1) The binding po- 
tential of magnetic impurity is strong enough to generate 
significant weight of the carriers of one spin at the im- 
purity site. This results in a large splitting between the 
majority and minority spin bands, and leads to quali- 
tatively different band structures from continuum model 
predictions. (2) The spin polarization of itinerant carriers 
increases with the interaction with magnetic impurity, J, 
for fixed doping concentration, x, and also increases with 
x for a fixed J. The intrinsic spin-orbit coupling effects 
reduce the degree of spin polarization and have relatively 
bigger effects on spin polarization at low x. The ability of 
manipulate the spin polarization by doping is attractive 
for magnetotransport applications. (3) For low doping, 
spin orbit and exchange couplings lift the degeneracies 
of the impurity bands resulting in an insulator with the 
chemical potential in the band gap. 

We describe DMS following the sp-d exchange 
model |15| | containing the Hamiltonian of the host semi- 
conductor in the clean limit and the effective exchange 
coupling between the localized and itinerant spins. Since 
the relevant energy scale is only order of a few hundred 
meV, it is of crucial importance to use a host semicon- 
ductor Hamiltonian which is precise close to the Brillouin 
zone center. To this end, we use the Effective Bond Or- 
bital Model (EBOM) developed by Chang EBOM 
has previously been applied to DMS within a contin- 
uum mean- field approximation ,251]. The difference in 
our model is the explicit inclusion of the localized char- 



FIG. 1: The impurity binding energy as a function of ex- 
change coupling constant J p -d- The experimental value is 
indicated by the dashed line. The inset shows the orbital- 
resolved local density of states at the impurity site with the 
fitted exchange coupling of J p _cj=-2.48 eV. The density of 
states are vertically shifted for clarity. 

acter of the exchange interaction. Our basis set is com- 
posed of the bulk state solutions at k=0, which arc Bloch 
sum of effective bond orbitals located at fee lattice sites. 
This basis set has the nice feature that the second order 
expansion of the matrix elements about the zone cen- 
ter yields the results of the k ■ p model. In addition, it 
reproduces the experimental band structure over wider 
ranges of the Brillouin zone. The impurities are placed 
at the fee lattice sites and the matrix elements are cal- 
culated in the effective bond basis. The applicability of 
this scheme in calculating impurity states has been well 
demonstrated 01 • The exchange couplings are responsi- 
ble for the hybridization of impurity states with itinerant 
electron states as well as the interaction between elec- 
trons at impurity sites. In general, it is possible that the 
contribution from electron-electron interaction can give 
rise to a dependency of exchange coupling constants on 
magnetic impurity concentrations. We believe, however, 
this dependency is negligible because of the large dielec- 
tric constant of GaAs and the long Fermi wavelength of 
typical DMS systems. The resulting Hamiltonian can be 
written as 

H — f 1 " 1 ^ r + E nm r* r 

^ ^ JnmSlCi n &nmCl7n ) (1) 

where i, and j are the fee site indices, / is the location 
of the substitutional impurities on the fee sites, <r is the 
Pauli spin matrix, and Sj is the local spin moment, which 
we treat as a classical object of magnitude 5=5/2. n and 
m are the combined indices of the spin and the bond or- 
bital, i.e., n — {t, a}, where in the current 8-band model 
r={|, 1} and a—{s,p x ,py,p z }. The first two terms on 
the right hand side of are the hopping and on-site 
potential of EBOM, respectively. The SO coupling is 
included in the onsite matrix elements E™ m . The last 
term in Eq. Q is the exchange interaction H ex , whose 
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FIG. 2: The valence band-edge splitting in the absence of 
spin-orbit coupling as a function of impurity concentration 
for various exchange couplings J p -d- The star symbols are 
the data from Ref. |23j|. The inset shows the valence band- 
edge splitting as a function of exchange coupling at a fixed 
impurity concentration of x = 0.037. 

coupling matrix J nm = —(n\H ex \m)(2/S) is to be de- 
termined as discussed below. Our formalism is general 
enough to include disorder Jljj; however in this paper 
we calculate the electron states for supercells made of fee 
lattices with a single magnetic impurity oriented in the 
z-direction (say). The magnetic impurity concentration, 
x, is defined as the inverse of the number of fee sites. 

The exact value of the exch ang e coupling is rather un- 
settled. Magneto-reflectivity |l9|. core-level photoemis- 
sion J2jj, resonant pho toemission spectroscopy |2l|, and 
magneto-transport [22] measurements found 2.5 eV, -1.0 
eV, -1.2 eV, -3.3 eV for the p-d exchange coupling, re- 
spectively. To circumvent this ambiguity, we determine 
the exchange coupling matrix by comparing the impurity 
bound state energy with experimental data. Because Mn 
atoms in GaAs are substitutional impurities with the site 
symmetry of the Ga sites [2(|, it can be readily shown 
that the exchange coupling does not mix different or- 
bital states. The exchange coupling of conduction elec- 
trons with Mn 3d electrons, J s -d = J T s,ts, are known 
to be very weak, generally of order of 0.1 eV 0. We 
found that our results are very insensitive to this cou- 
pling. Hereafter, we set the coupling to the conduction 
electrons to 0.1 eV. 

In Fig. ^ we show the binding energy of a single impu- 
rity as a function of the exchange coupling. The bound 
state energies were extrapolated from supercell calcula- 
tions of 13 3 to 15 3 fee lattice systems. With the best 
fitted exchange coupling J p -d = -2.48 eV, the impurity 
bound state energy agrees with experimental observa- 
tion [53 of 0.113 eV. The orbital-resolved local density 
of states, presented in the inset of Fig. ^ shows that 
the bound states are in almost definite spin states and 
split into orbital angular momentum substates of m 2 =- 
1,0, and 1. This agrees with recent calculation by Tang 
and Flatte based on atomic orbital tight-binding Green's 
function method It is worth noting that there is a 
critical magnitude of the exchange coupling ^1.8 eV, be- 
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FIG. 3: (a) The valence band structures of Gao.963Mno.o37As 
calculated from EBOM model with exchange coupling J p -d 
= -2.48 eV. The Fermi level is depicted by the horizontal 
line. The spin-resolved density of states is drawn on the left 
(right) for majority (minority) spins, (b) Same figure as (a) 
except calculated in the absence of SO coupling, (c) kp band 
structure with mean- field exchange coupling (3No = -1.2 eV 
and x = 0.037. 



low which the valence holes are not bound to impurities. 

The applicability of our model to dense impurity con- 
centration regime is supported by a comparison with ex- 
isting ab initio DFT calculations [HHj. In the DFT 
calculation by Sanvito et al. [23j, the band-edge split- 
ting of Gai-^Mn^As was calculated as a function of the 
impurity concentration. We performed equivalent super- 
cell calculations using our model Hamiltonian, turning 
off the SO coupling for comparison purposes. In Fig. [21 
we present the band-edge splitting of the majority and 



minority valence bands. At J p 



J-, 



= -2.40 cV, 



our calculation shows an excellent agreement with the ex- 
isting DFT calculation in the region x ~ (0.01, 0.04) p3| . 
The fitted exchange coupling differs only by 0.08 eV from 
the fitting to the single impurity bound states. This is 
an encouraging result since it confirms that our model 
is applicable in a wide range of impurity concentrations. 
The calculated band structure with the fitted coupling 
parameter, Fig. Gib), also agrees very well with previous 
DFT calculations j2J, [24| . The nonlinear dependency of 
valence band-edge splitting on the exchange coupling, in- 
set of Fig. [21 shows a clear difference between our model 
and the virtual crystal approximation ( VCA) , which pre- 
dicts a linear relation when SO coupling is not present. 

Having established the model, we turn to assessing 
the effects of SO coupling on the electronic structure. 
Most of DFT calculations have predicted half-metallicity 
in Gai-zMn^As [23l Izif and this prediction has drawn 
much research interest in regards to possible application 
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FIG. 4: Density of states Gai-^Mn^As. (a) Magnetic impu- 
rity concentration x — 0.037. Inset shows the polarization as 
a function of exchange coupling, (b) Exchange coupling, J p -d 
= -2.0 eV. Inset shows the polarization as a function of im- 
purity concentration. Circle, square, and triangle correspond 
to J p _<j=-1.5 eV, J p _d=-2.0 eV, and J p -d~-2A8 eV, respec- 
tively. For both (a) and (b) the upper panel is the majority 
spin density of states and the lower panel is the minority spin 
density of states. The energy is relative to the Fermi energy 
located at zero. 



to the spin devices |2|. When the spin degree of free- 
dom is separated from orbital degrees of freedom, it is 
natural to have spin-split bands in the presence of mag- 
netic background. In Fig- El (b) we show the band struc- 
ture of Gao.963Mno.o37As when SO coupling is not con- 
sidered. The separation of orbitally degenerate majority 
and minority spin bands agrees very well with DFT cal- 
culations |23| . On the contrary, in k • p approximations, 
where the impurities are treated within a virtual crystal 
approximation, the exchange interaction is scaled by im- 
purity concentration and ends up introducing a relatively 
weak potential. Fig. [31(c) shows a typical k-p bands with 
mean-field exchange coupling f3N = -1.2 eV The ex- 
change coupling splits the heavy- and light-hole bands as 
little as 0.1 eV and the mixing of light- and heavy-hole 
bands leads to weak spin polarization. Our calculation 
resolves the conflict between above two predictions. Tak- 
ing into account both SO coupling and local exchange in- 
teraction, we found that electron states of Gai-^Mn^As 
with dense impurity concentration, i.e. x > 0.01, is in- 
deed half-metallic. Fig. 0(a) illustrates that the splitting 
due to exchange coupling is much larger than the SO cou- 
pling and the low states are fully spin polarized. 

To explore the dependency of the spin-polarization, 
we present the spin-resolved density of states for vari- 
ous exchange couplings and impurity concentrations in 
Fig. 0] When exchange coupling is not strong enough 
to induce impurity bound states, e.g. J p -d ~ -1-5 eV, 
the minority spin has substantial contribution to the po- 
larization. As the exchange coupling gets stronger the 
polarization becomes bigger, reaching near 100% at J p ~d 
~ -2.48 eV. The dependency of spin polarization on the 
impurity concentration is of more practical interest since 
it can be controlled in sample growing processes. The 
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FIG. 5: The valence band structures of Gao.99sMno.002 As. 
(a) without spin-orbit SO coupling, (b) with SO coupling. 
The dotted line shows the location of the Fermi level. Note 
that in (b) the Fermi level lies in the band gap generated 
largely by SO effects. 

increase in the impurity concentration contributes to the 
spin-polarization by enhancing the splitting between the 
majority and minority bands. We found the spin po- 
larization monotonically increases within the experimen- 
tally accessible impurity concentration range x < 0.1. 

At low impurity concentrations, the SO coupling can 
be comparable to the width of the impurity bands and 
this can result in a novel type of insulating phase. The 
mechanism by which such an insulator is formed is dis- 
cussed below. A single orbital is split into two spin polar- 
ized bands by the exchange interaction. In the absence 
of charge compensation the majority spin band is com- 
pletely filled while the minority band is empty. Since 
the chemical potential lies in the middle of the spin-split 
bands, the system is insulating. In a real semiconductor, 
on the other hand, with several degenerate orbitals in 



the valence band, the degeneracy can be lifted by two ef- 
fects arising from the exchange coupling as well as from 
the spin-orbit coupling. The exchange coupling acting 
alone splits the orbitals into partially occupied majority 
bands and empty minority bands from which one would 
have concluded that the system is metallic as seen in 
Fig-EJa). However, the presence of spin-orbit coupling 
further splits the bands and opens up a gap. For suf- 
ficiently narrow impurity bands, the chemical potential 
then lies within the gap and results in an insulating state 
as seen in Fig.[SJb). In strongly SO coupled systems, such 
as InSb, one might be able to find the insulating phase 
driven by the SO coupling even with moderate disorder 
effects. Disorder not only broadens the impurity bands 
but also reduces the band splitting, resulting in the merge 
of impurity bands into the main valence bands in case of 
Gai^Mn^As |l8|. 

In summary we presented a tight-binding model which 
properly includes SO coupling effects in Gai^Mn^As. 
Our model reproduces the experimentally observed im- 
purity binding energy in dilute impurity limit and agrees 
well with first principle calculations in dense impurity 
limit. We found the impurity bound states are split due 
to the SO coupling. In dense doping regime, the spin- 
polarization increases with magnetic impurity concentra- 
tion and almost fully spin-polarized metallic states can 
be achieved. We observe impurity bands formation at 
low impurity concentrations and predict that strong SO 
coupling can induce a novel insulating phase. 
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